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Kshirsagar-Tan independence property of 
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A new independence property of univariate beta distributions, related to the results of Kshirsagar 
and Tan for beta matrices, is presented. Conversely, a characterization of univariate beta laws 
through this independence property is proved. A related characterization of a family of 2 x 2 
random matrices including beta matrices is also obtained. The main technical challenge was a 
problem involving the solution of a related functional equation. 
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1. Introduction 

Characterizations of probability distributions by transformations preserving indepen- 
dence have been investigated by many authors. Celebrated theorems of Bernstein or 
Darmois and Skitovitch (for the normal law), of Lukacs (for the gamma law) and of Fisz 
(for the exponential law) are of this form. For beta laws, characterizations of this na- 
ture are typically related to neutrality properties of the Dirichlet distribution; see Fabius 
(1973), James and Mossiman (1980), Geiger and Heckerman (1997) and Bobecka and 
Wesoiowski (2007). Matrix variate analogues of such characterizations are often much 
more demanding; see, for instance, characterizations of Wishart matrices due to Olkin 
and Rubin (1964), Casalis and Letac (1996), Geiger and Heckerman (2002), Bobecka 
and Wesoiowski (2002) or characterizations of Wishart and matrix variate generalized 
inverse Gaussian matrices obtained in Letac and Wesoiowski (2000) or Massam and 
Wesoiowski (2006). 

Our general goal is to develop characterizations related to independence properties 
of matrix variate beta distributions. A characterization of this type, related to matrix 
versions of neutrality properties, has recently been given by Hassairi and Regaig (2006). 
They extended to matrix random variables the result for univariate beta laws obtained 
in Seshadri and Wesoiowski (2003). Here, we consider independence properties of beta 
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matrices discovered by Kshirsagar (1961, 1972) and Tan (1969), which are essentially dif- 
ferent than neutrality. This is carefully explained in Section 2. The study of the easiest 
possible case of 2 x 2 matrices led us to a new independence property for three univariate 
independent beta variables (for a recent review of the univariate beta distribution, one 
can consult Gupta and Nadarajah (2004)). The transformation preserving independence, 
we propose, seems to have no direct connection to neutrality (see the beginning of Section 
3 and Section 5). The main results of this paper, given in Section 3, are characterizations: 
(1) of univariate beta laws through this new independence property; (2) of a family of 
2x2 matrix variate distributions satisfying Kshirsagar-Tan independence properties. 
Somewhat unexpectedly, this latter family happens to be wider than the class of 2 x 2 
matrix beta distributions. The proofs are based on a solution of a related functional 
equation. This equation is solved, under technical smoothness conditions, in Section 4. 
Concluding, in Section 5, we discuss similar characterizations for Wishart-type distribu- 
tions and relations of our results to neutrality and to a representation of the Dirichlet 
distribution, as well as to uniqueness of solutions of some stochastic equations. 



2. Independencies for beta matrices and beta 
distributions 

Let Vn denote the Euclidean space of n x n symmetric matrices (with real entries) with 
the inner product defined as a trace of the product. The Lebesgue measure on V„ is 
fixed by assigning the unit mass to the unit cube in this space. Let C V„ denote 
the cone of positive definite symmetric real n x n matrices. Additionally, denote by 
M n^m the space of tt. x to real matrices. Let e„ be the n x n identity matrix. Also, let 
2?„ = {x G V+ : e„ — cc G V+}, an analogue of the interval (0, 1) in V+. 

The matrix variate beta distribution /3„ {p, q) on I?„ is defined by the density (with 
respect to the Lebesgue measure on Vn) of the form 

, (detu)P-("+i)/2(det(e-w))«-("+i)/2 

fm = —. — r iv„{u), 

B„{p,q) 

where p,q> ^^^^ and _B„ (p, q) is the n-dimensional Euler beta function defined in terms 
of 71-dimensional Euler gamma functions as 

o , ^ r„(p)r„(g) 

Bn[p,q) = -fT-^ — 
Recall that (see, e.g., Muirhead (1982), page 62) 

i-l 



r„(p) = 7r"("-i)/4[]r 



Let X be a beta n x n random matrix /3„(p, q). Let T = [tij] be an upper triangular 
random matrix defined by TT' = X. Kshirsagar (1961, 1972) (see also Muirhead (1982), 
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Chapter 3.3 or Gupta and Nagar (2000), Chapter 5) proved that i = are 
independent beta variables, i| ^ Bi{p — g), i = 1, . . . , n. In the course of the proof, he 
considered the block partitioning of T with the dimensions of the blocks lxl,lx(n— 1), 
(n — 1) X 1 and (n — 1) x — 1), as 



T- 



tii t' 
T22 



and showed that 



tii,T22 and v = 



(e„-i — T22T22) 



-1/2, 



are independent. Moreover, it easily follows from his proof that ^22^22 has the matrix 
beta distribution /3„_i(p— i, q) and v has the density /v(x) oc (1 — x'x)'?~("+i)/2/^^_^(x), 
where 5„_2 is a unit ball in R"^i. This result was extended in Tan (1969), who considered 
the block partitioning of the matrix X ^ /3„(p, q) with the dimensions of the blocks r x r, 
r X s, s X r, s X s, 



X 



Xu X12 
X21 X22 



and proved that 



^11, -'^2-1 = -'^22 - -'S^21-'^ll^-'^12, V^{es-X2.l) ^^'^X2l{Xij'-{er - Xii) ^)^^^ 

arc independent random matrices. Moreover, it appears that Xn is a beta matrix (3r{p, q), 
X2-1 is a beta matrix /3s (p— §,?) and V has the density 



fvix) cx [det(e, - x'xW-^''+^^/^Ih,,Ax), 



where Tir 



{x e Ms,r ■ er - x'x G V,,t}. 



Our aim is to study a converse problem, that is, we want to know if the independence 
conditions mentioned above characterize the matrix variate beta distribution. Since, in 
the general case, the problem seems to be very hard, here we study only the case of 2 x 2 
matrices. Let 



X = 



Xii X12 
X12 X22 



be a beta matrix (32{p,q)- It then follows from the results of Kshirsagar (1961) and Tan 
(1969) that 



-'^11, -^^2-1 = ^22 7; — ,Vi = 

All 

A12 



Xi2 



■\/(l — -^ll)(-'^ll — -^^11-^22 + -^^12) 
-'^12 



-^^22, -^^1-2 — -'^ll 77 — , V2 — — , 

^22 v/(l-^22)(^22- ^11^22 



xh) 



are independent, (2.1) 
are independent (2.2) 
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and both the triplets have the same distribution: Xu and X22 have the beta distribution 
Bi{p,q), ^21 and X1.2 have the beta distribution Bj{p— ^,9), Vi and V2 have the 
beta distribution Bj{^,q— i). It appears that the independence conditions (2.1) and 
(2.2) characterize a family of distributions which is wider than that of the beta matrix 
distributions. 

First, instead of (2.1) and (2.2), we consider slightly weaker conditions, in the sense 
that the third clement of both the triplets is squared. This allows us to rephrase the 
independence property completely in terms of univariate beta variables. As a matter of 
fact, an even more general property holds. This will be explained now. 

Let 

H = {{xi,X2,X3) e (0,00)^ :min{xiX2, (1 - a;i)(l - X2)} > X3}. 
Denote by B{p,q,r) the probability distribution on H with the density 

f{xi,X2,X3) (X {xiX2 - X3f^'^[{l - .Ti)(l - X2) - Xa]'^^'^ x^^^^'^ I H (xi , X2, X3) , 

where p, q, r are positive numbers. 
On H, define two mappings 

( X\X2 ^ X3 Xj, \ 

'\l}^{xi,X2-,x-i)^\Xi, ,- ■ , z = l,2. (2.3) 

V Xi (\-Xi)(Xi-X\X2+X-i) } 

Note that "0^, 1 = 1,2, are bijections of H onto (0, 1)^. Moreover, it can be checked that 
^' = V'2 o V'r^ ■ (0' 1)^ ~^ (0' 1)^ is involution (i.e., ^' = 5""^) and for any (yi,y2, ys) € 

(o,l)^ 

*(yi,y2,2;3) 

= ^2/2 + (i-yi)(i -2/2)^3, 

yiy2 

y2 + (1 - yi)(i - 2/2)2/3 

Let (Xi,X2,X3) be a random vector having the distribution B{p,q,r). A standard 
calculation using the Jacobian shows that Tjji{Xi, X2, X3) is a random vector with inde- 
pendent beta components 

^Pi{Xi,X2,X3) - Blip + r,q + r)® Bi{p, q + r) ® Bi{r, q), i^l,2. 

This observation can be rewritten as an independence property of a transformation of 
invariance beta variables. 

Theorem 1. // (Yi, ^2, 13) ~ Bi{p + r,q + r)® Bi{p, q + r)® Bi{r, q), then 



(2.4) 



2/12/3 



[1 - (1 - yi)y3][y2 + 2/3(1 -^2)] 



i"2 + (l-n)(l- 1^2)^3, 
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>2 + (1 - n)(l - Y2)Y3 ' [1 - (1 - Y,)Y3][Y2 + 1:3(1 - Y2)] 

3. Characterizations 

First, we study a characterization of beta laws which is a converse of Theorem 1. As 
aheady mentioned in the Introduction, to the best of our knowledge, all characterizations 
of beta variables existing in the literature which are based on transformations preserving 
independence can be related to neutrality properties of the Dirichlet distribution. The 
characterization of beta variables we give below is based on the independence of the 
components of the random vector defined in Theorem 1 . It seems to have no connection 
to analogous results based on neutralities. This issue will be explained more carefully in 
Section 5. 

Theorem 2. Let Y — (Yi, 12,13) be a (0,1)'^ -valued random vector with independent 
components having strictly positive, continuously differentiable density on (0, 1)'^. Let ipi 
and ip2 be the transformations defined in (2.3) and = -02 o ip^^ . Let 

z = (Zi,Z2,Z3) = *(ri,r2,i"3) 

also have independent components. 

There then exist positive numbers p, q, r .such that 

Y = Z^Bi{p + r,q + r)<g)Bi{p,q + r)®Bi{r,q). 

Proof. Note that Y = -01 (X) and Z =-02 (X) , where X is a random vector assuming 
values in H. The Jacobians oi ipi, i = 1,2, are, respectively, 

1^.1 = 7^ TT^ — ^' ^=1'2- 

(1 - Xi)(Xt - X1X2 + xa) 

Thus, the density / of X can be expressed as 
/(x) = ^ 



(1 -xi){xi -X1X2 +X3) 
2 - 3^3 

Xl J " ^ \{\ - xi){xi - X1X2 + x^) 



r t -.r I X1X2 X3 



(3.1) 



(l-.T2)(x2-xia;2+a;3) 



(1 - a;2)(a;2 - X1X2 + X3) 



(3.2) 
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We multiply (3.1) by 2:3 and define 

gi(.T) =ln(/yi(j:)), g2{x) ^Inifv^ix)), g3(x) ln(a-/y3 (a;)), 
g4(.T) =ln(/zi(a;)), 55(0;) = ln(/z2 (x)), ge{x) ^Inix/z^ix)). 
We then obtain 

51(2/1) +52(2/2) +53(2/3) =.94(21) +.95(^:2) +.96(23), (3.3) 
where ipiix) = y and -02 (x) = z. Note that y e (0, 1)"^ and that 

21 =2/2 + (l-2/i)(l- 2/2)2/3, 22 



23 = 



2/2 + (1- 2/0(1-2/2)2/3' 
2/12/3 



(1 - (1 - 1/1)2/3) (2/2 + 2/3 - 2/22/3) ' 



Note that equation (3.3) is the one we solve in Proposition 1 (Section 4). Hence, we 
conclude that there exist constants p, q, r such that 

JyA^) = fzA^) = BixP+^-'C^ - x)''+--\ 
fY,{x) = fzA^)^B2xP-\l ~ xy+'''-\ 
fYA^)^fzs{^)=B3x''~Al~x)'^-\ 

where the Bj's are some normalizing constants. The integrability condition implies that 
p, q, r are positive. □ 

An alternative formulation of the above theorem is as follows. 

Theorem 3. Let X= (Xi, X2, X^) be an H -valued random vector having strictly pos- 
itive, continuously differentiable density. //■i/;i(X) and -02 (X) have independent compo- 
nents, then there exist positive numbers p,q,r such that ^B{p,q,r) . 

If, instead, we use original independencies which hold for the beta 2x2 matrices, that 
is, (2.1) and (2.2), the functional equation of Proposition 1 (Section 4) leads to a family 
of distributions on P2 which is wider than that of beta matrix distributions. 



Theorem 4. Let X be a 'D2-valued random matrix having a density which is continuously 
differentiable and strictly positive on I?2- Assume that (2.1) and (2.2) hold. There then 
exist constants a,b,c> such that the density of X is of the form 
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Proof. Using (2.1) and (2.2) similarly as in the proof above, we get for the density fx 
the following two representations: 

1 



a;ii)(xii - 2:iia;22 + xj^) 

„2 



xii J V v/(l - xii)(xii - xiia;22 +0:^2) 
1 



- X22)(x22 - a;iia;22 + a;?2) 



X11X22 - a;f2 \ J. / a;i2 



x/4(x22)/5^^^^ ^ /O, ^ 

V 2:22 / Vv(l-a;22)(x22-a;iia;22+a;T2) 

where fi, i = 1,2, ... ,6, are respectively densities of Xii,X2.i, Vl,-''^22,-'^^i-2, ^2- We now 
proceed similarly as in the proof of Theorem 2, separately in two cases: X12 > and 
Xi2 <0. 

For X12 > 0, denote \fu — xyi- The above equation then takes on the form 
UxM^-^^^^\h 



\ Xxx ) \{l - Xx\){xxx ~ XXXX22 ^ u) 
^/4(x22)/5f^^ii^Y7a^ 



V 2:22 /' V (1 - 2:22) (2:22 - a;iiX22 +m) 

where /i(a;) = \fxji{^/x), i = 3,6. Thus, from Proposition 1, it follows that there exist 
positive (due to integrability) constants a,h,c such that for x g (0, 1), 

h{x) = /4(x) = Bi.t'^+^-1(1 - x)''+^-i, 
h{x) = ]^{x) - B22-''-i(l - a;)^+'=-i, 
/3(x) = /6(x)- 53x2^-1(1 -x^)"-!. 

For a:;i2 < 0, denote \fu~ —x\2. The above equation then takes on the form 
A(xn)/2f^^ii^)/3 



V Xx\ ) \{^ - Xx\){xxx ~ Xx\X22 + u) 
, I w / Xx^X22 i ( U 



V 2:22 / V (1 ~ 2:22) (3:22 - •i;ilX22 + It), 

where fi{x) = \/xfi{—y/x), i = 3,6. From Proposition 1, it follows that for x G (—1,0), 
h{x) = h{x) = i?3(-x)2=-i(l - x^ f-\ 
Thus, /3(2:) = /6(x)=i33k|2--i(l-a;2)''-i/(_i,i)(x). 

Finally, the result follows from the representation of fx in terms of /i,/2,/3. □ 
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Remark 1. Note that the family of matrix variate distributions defined by (3.4) includes 
the beta matrix distribution: (3.4) for c = ^ is the density of the beta matrix distribution 
/?2(a+i,6+i). 

4. Functional equation 

In this section, we solve the functional equation which was essential for proving the 
characterizations derived in Section 3. 

Proposition 1. Let gi, i = 1,2,3, be continuously dijjerentiable functions on (0,1) sat- 
isfying 

51 (yi) +32(^2) +53(^3) 

= 54(2/2 + (1 - yi){l - 2/2)2/3) + 55 (^TT^^, r-) (4-1) 

V2/2 + (l -2/i) (1-2/2)2/3/ 
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2/12/3 



(1 - (1 - 2/i)2/3)(2/2 + (1 - 2/2)2/3) 



for any yi,?/2,2/3 S (0, 1). 

There then exist constants a, f3,"f, Ai, i = 1, . . . , 6, such that Ai + A2 + A3 = Ai + A5 + 
Aq and 

g,(.T) = (a + 7)ln.T + (/3 + 7)ln(l-x) + Aj, i = l,4, (4.2) 
gi{x) = a\nx + {f3 + -f)\n{l - x) + Ai, i = 2,5, 

gi{x) = jinx + p\-a{l — x) + Ai, i = 3,6. 

Proof. Let s = 5(2/1,2/3) = 1- (l-2/i)2/3, * = ^(2/2,2/3) = 2/2 + (1-2/2)2/3- Taking derivatives 
of (4.1), once with respect to j/i, once with respect to 1/2 and once with respect to j/s, 
yields the following three equations: 

g'livi) = -(1 - 2/2)2/32/4(1 - (1 - 2/2)5) (4.3) 

V2t I f 2/12/2 \ 2/3(1-2/3) , / 2/12/3 \ 

+ [1-(1- 2/2)5]2^^ 1^1 -(1-2/2)5/ ,2t .St 

52(2/2) = 5.9^(1 -(1-2/2)5) (4.4) 

2/1(1 - 5) , / 2A2/2 \ _ 2/1(1-2/3)2/3 , / 2/12/3 ' 

^ [1-(1- 2/2)5]2^^W- (1-2/2)5/ 5t2 5i 

53(2/3) = (1 - 2/i)(l - 2/2)54(1 - (1 - 2/2)s) (4.5) 
2/1(1 -2/1)2/2(1-2/2) , / 2/12/2 



[l~(l-2/2)5]2 ■95Vl_(l_y,), 
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, 2/i[y2 + (1-^2)2/3(1-5)] ifyiy^" 



Denote 

Gjix)=xg'^{x), (0,1), j = l, 2,3,4,5,6. 

Then, upon multiplying (4.3), (4.4) and (4.5), respectively, by yi, 1/2 and 2/3, we get 

t yiy2 ^ ^ ^~^^ Co^ ^^^^ 



i-(i-2;2)s yi- (1-2/2)5/ s \ si 

G2(y2) - , , G4(l - (1 - 2/2)5) (4.7) 

1- (1-2/2)5 

_l_ (1-2/1)2/3 / 2/12/2 \ _ 2/2(1-2/3) ^^ 2/12/3 



1- (1-2/2)5 °Vi- (1-2/2)5; t °V 5t 



and 



G3(2/3) = ^}^|^^G4(1- (1-2/2)5) 

(1 -2/2)(l - 5)^ / 2/12/2 \ ^ 2/2 + (1 - 2/2)2/3(1 - ^) (^Y 



(4.8) 

i-(i-2/2)5 ""vi- (1-2/2)5; ' Si ""V St ;■ 

Letting j/3 ^ in (4.6), we eoncludc that the limit 

hmGaf^ 

y.-s^o \ st 

exists. We denote it by 6*6(0). (4.6) then yields 

Gi(2/i) = G5(2/i) + G6(0). (4.9) 

Similarly, (4.6) for 2/3^0 gives 

62(2/2) = G4(2/2)-G6(0). (4.10) 
Multiplying (4.7) by 1 — j/2 and adding the resulting equation to (4.8), we get 

(1 - 2/2)62(2/2) + 63(2/3) ^ (1 - 2/2)64(1 - (1 - 2/2)5) + ^^l^l^y'^y^ G, . (4.11) 

Letting j/i 1 in (4.11), we obtain 



(1 - 2/2)62(2/2) + 63(2/3) = (1 - 2/2)64(2/2) + 2/266 
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Thus, by (4.10), we have 

-(1-^2)06(0) + 63(^3) =2/2^6 I 
and letting j/2 — > 1 , we get 



G3(y3) = G6(y3)- (4.12) 

From (4.11), we have 

^ f ym\ _ 1- (1-2/1)2/3 



st J 1 — (1 — 1/2)3 

(4.13) 

X [(1 - 2/2)62(2/2) + 63(2/3) - (1 - 2/2)64(1 - (1 - 2/2)s)]. 

We now plug (4.13) mto (4.6) and (4.7). Multiplying the resulting equations by (1 — 
yi)j/3 and —t, respectively, and then adding them, we get 

(1 - 2/1)2/361(2/1) - 62(2/2) - (1 - 2/3)63(2/3) = -564(1 - (1 - 2/2)5). (4.14) 

Letting 1/2 0, we obtain 

(1 - 2/1)2/361(2/1) - 62(0) - (1 - 2/3)63(2/3) = -.s64((l - 2/1)2/3). (4.15) 

Note that the hmit linij^^o 62(2/) exists. We denote it by 62(0). Thus, (4.10) implies that 
the limit limy^o64(y) also exists. We denote it by 64(0). (4.10) then implies 64(0) = 
62(0) + 65(0). Subtracting (4.15) from (4.14) and again using (4.10), we get 

64(0) - 64(2/2) - s64((l - 2/1)2/3) = -564(1 - (1 - 2/2)5). 

Let X = (1 — 2/1)2/3 G (0, 1). Then, for any x,y — 1/2 (z (0, 1) from the above equation, we 
obtain 

64(0) - 64(2/) - (1 - a;)64(x) ^ -(1 - x)Giiy + V^)- (4.16) 
Changing x into y and y into x, we get 

64(0) - Gi{x) - (1 - y)Gi{y) = ~{l ~ y)Gi{y + x~ yx). (4.17) 

Multiplying (4.16) by 1 — 2; and (4.17) by 1 — x and then subtracting the resulting 
equations, we get 

-2/64(0) - (1 - 2/2)64(2/) - (1 - 2/)(l - x)Gi{x) 
= -.t64(0) - (1 - x)Gi{x) - (1 - x){l - y)Gi{y)- 
Hence, for any x,y<^ (0, 1), 

--64(0) + ^64(2/) = --64(0) + i^64(x) 

y y XX 
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and by separation of variables, we get 

--G4(0) + i^G4(x)=C, 

X X 

where C is a constant. Thus, for any x € (0, 1), 

G4(.t) = -C+[C7 + G4(0)]-^ (4.18) 

1 — X 

and we get 

54(2;) =G4(0)lnx- (G + G4(0))ln(l -x) +v44- 
From (4.18) and (4.10), it follows that 

G2(x) = -(G + G6(0)) + [G + G4(0)]-^ (4.19) 

\ — X 

and we obtain 

gi{x) = (G4(0) - G6(0)) Inx - (G + G4(0)) ln(l ~ x) + A2. 
Plugging (4.18) and (4.19) into (4.15), wc get 

(1 - 2;i)2/3Gi(yi) + G6(0) - (1 - yz)G^[y3) = -(1 - Vi)yzC. 
Hence, for any yi,?;3 e (0, 1), 

(1 - 2/i)Gi(yi) + (1 - 2/i)G = + i^G3(2/3) 

2/3 ys 

and again using separation of variables, we have 

{l^y)G,{y) + {l-y)C = D 

and 

^6(0) ^l-V, 

y y 

where D is a constant. Thus, 



-G3(2/)=A 



1-2/ 

and 



Gi(y) = -G+-^ (4.20) 



G^{y) = -D+[D + G^m-^. (4.21) 

1-y 



Hence, 

gi{x) =^ (-G + £>) Inx - i:'ln(l - x) + Ai 
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and 

53(2;) = GeiO) \nx~iD + GeiO)) ln(l ^ x) + A3. 
Using (4.9), we obtain 

G5{y)^-{C + GeiO)) + -^ (4.22) 

and using (4.12), we get 

Ge{y)^-D+[D + Gem-^. (4.23) 

Hence, 

g5{x) = i-C + D - G6{0))liix - Dlnil - x) + A5 

and 

ge{x) = Ge{Q) \nx ~ [D + Gq{Q)) ln(l - x) + Ag. 
Inserting the formulas obtained for g^, i = 1, . . . , 6, into (4.1), we get that it must be 

GA{i)) = D-C. 

Thus, we get (4.2) where 

a = D-G-GQ{Q), P = -D-Ge{0), -f = Ge{0). □ 

5. Concluding remarks 

The transformation ^! preserving a product beta distribution considered in this paper 
is defined in the thrcc-variate situation. Describing its higher-dimensional version is an 
open problem. 

On the other hand, it is natural to seek a multidimensional analogue of Theorem 4 by 
referring to matrix variate versions of independencies, as described in Section 2, for all 
or several block partitionings of the original matrix. 

Recall that characterizations of Wishart matrices using independencies of {Xn, X12) 
and X2-1 for all or several block partitionings of the n x n random matrix X with n> 2 
were proven in Gcigcr and Hcckcrmann (2002) and in Massam and Wesolowski (2006). 
It is interesting to note that for 2x2 matrices, similarly as in Theorem 4 here, these 
independencies characterized wider families of random matrix distributions; see Geiger 
and Heckcrman (1998) or Letac and Massam (2001). 

Also, note that if X is Wishart with a diagonal matrix parameter, then for any block 
partitioning of X , the following independence properties hold: Xn, X2iXfj^Xi2 and X2.1 
are independent and X22, ^12^22^^21 and X1.2 arc independent (sec, e.g.. Example 3.14 
in Muirhead (1982)). Recently, a characterization using these independence properties 
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for 2 X 2 matrices was proven in Seshadri and Wesolowski (2007). Analogously to our 
Theorem 4, it was shown there that the two independence conditions characterize a 
family of random matrices which is wider than that of Wishart matrices with diagonal 
matrix parameter. Again, the problem in higher dimensions remains open. 

It is natural to compare the transformation ^' considered here to neutrality properties 
of the Dirichlet distribution (being consequences of the gamma variables representation), 
which define transformations preserving independence of beta variables. In particular, 
the one defined through the concept of complete neutrality (see Connor and Mosimann 
(1969), James and Mosimann (1980)) is worth analyzing. In the three- variate case which 
we are interested in here, the respective characterization has the following form. Let 
Y = (Yi, I2, ^3) be a (0, l)^-valued random vector with independent components. Let 

Z=f 7 , (l-yi)y2 — (^_Y)n^Y2)Y3] (5.1) 

also have independent components. There then exist positive numbers p,q,r,s such that 
Y Bi{p,q + r + s) <E) Bi{q, r + s) <E) Bi{r, s) 

and 

Z - Bi{p, s) (g) Bi{q,p+s) (g) Bi{r,p + q + s). 

Let us emphasize that the transformation defined by (5.1) and ^ (sec (2.4)) used in 
Theorems 1 and 2 are essentially different, though they may seem similar at first glance. 

We conclude this section with two diverse and somewhat unexpected consequences of 
Theorem 1, one related to a representation of the bivariate Dirichlet distribution and the 
other related to stochastic equations and perpetuities. 



Remark 1. Consider a random vector (VFi, W2) such that 

Wi,Y:-^]-Bj{p,q + r)®Bj{r,q). 

Equivalcntly, (Wi,W2) has the Dirichlet distribution D{p,r,q). By Theorem 1. looking 
at the last two coordinates of the random vector ^'(11,12,^3), wc obtain the following 
representation: 



with 



(VFi,W2) = C/(Fi,l/2) 



C/=____^4___, v,=Y2 and ¥2= ^^'^'^^^ 



Y2 + {1-Yi){l-Y2)Y3' l-{l-Yi)Y3 



It can easily be seen that ( Vi + V2 ) [/ ^ Bi{p + r,q), Vi ^ Bj (p, q + r) and V2 ^ Bjj (r, p + 
q + r) are independent, where Bij{r,p + q + r) denotes the second type of beta distribution 
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defined by the density 

^'(^) = Bir,p + q + r){l + x)P+i+^- 

Also, note that {^^^,UVi, j^^) = ^(Yi, F2, ili)- Consequently, by Theorem 1, we 
conclude that '^i^y^ ^ Bi{p + r,q + r) and U{Vi,V2) ^ D{p,r,q) are independent. 

Remark 2. The results of Theorem 1 or 2 can be reviewed in the setting of stochastic 
equations. Namely, by comparing the respective coordinates of Y and \1/(Y), we observe 

that 

I. The stochastic equation (for unknown R) 

R^AR + B, (5.2) 

where {A,B) = (—(1 — ^2)^3,^2 + (1 — ^2)^3) and R are independent, has a solution 
R = Yi; 

II. The stochastic equation (for unknown S) 

S = CS + D, (5.3) 

where (C, D) = i^"^—^^^, l^JA^Xiy^"^ and S are independent, has a solution 5 = 

III. The stochastic equation (for unknown T) 

T^aT + b+^, (5.4) 

where (a,6, c) = (^, "'"^^'^^^^'"'^^ , — ^^y^ ^^^ ) and T are independent, has a solution 

It follows from the theory of perpetuities (see, e.g., Vervaat (1979) and Goldie and 
Griibel (1996)) that for equations (5.2) and (5.3), these solutions are unique. However, it 

is not known if T = ^ is the unique solution of (5.4). Possibly, the reference most relevant 
to this uniqueness problem is the paper by Chamayou and Letac (1991), where many 
examples of stationary distributions for compositions of random functions are considered. 
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